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Abstract: We study the possibility of evading astrophysical bounds on light 
pseudoscalars. We argue that the solar bounds can be evaded if we have a 
sufficiently strong self coupling of the pseudoscalars. The required couplings 
do not conflict with any known experimental bounds. We show that it is 
possible to hud a coupling range such that the results of the recent PVLAS 
experiment are not in conflict with any astrophysical bounds. 


1 Introduction 

In a recent paper, Zavattini et ah P] have reported a rotation of polarization 
of light in vacuum in the presence of a transverse magnetic held. If we 
interpret this rotation in terms of the coupling of a light pseudoscalar particle 
to photons, 

Un = > ( 1 ) 

we hnd that the allowed range of parameters is 1 x 10® GeV < < 

6 X 10® GeV and 0.7 meV < < 2 meV [T], where is the mass of the 

pseudoscalar. This range of parameters does not conhict with any laboratory 
bounds [2]. However these values are ruled out by astrophysical considera¬ 
tions jS] if we assume that the pseudoscalar is an axion. For a review of 
axion physics see, for example. Ref. jlj. The most stringent astrophysical 
limit comes from SN1987A which suggests that the mass of the axion cannot 
be greater than 0.01 eV [Sj. Assuming that this particle is an invisible axion, 
this in turn implies an upper bound on the coupling, < 2 X 10 ^^GeV \ 
upto a model dependent factor of order unity. The limit on the axion mass 
arises by demanding that axion emission should not lead to too much energy 
loss from the core. Similar but less stringent limits can be obtained by consid¬ 
ering energy loss from the core of the sun. The coupling, Eq. ^ also leads to 
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several interesting astrophysical polarization effects jB]. For standard axion, 
its mass and coupling to photons are both related to the Peccei-Quinn scale 
[7j. It is clearly of interest to see if the astrophysical bounds can somehow be 
evaded [S]. All the astrophysical bounds assume that the pseudoscalar cou¬ 
plings are so small that once produced it will freely escape from the source, 
which may be the sun or a red giant or a supernova. In the present paper 
we examine whether the pseudoscalar particles might be trapped inside the 
sun. 


2 Trapping Pseudoscalars 


We assume the following interaction lagrangian of pseudoscalars, 


Hi 




( 2 ) 


Here we have also included the self coupling of the pseudoscalars. The value 
of the pseudoscalar coupling to photons will be taken from data p. We do 
not assume 0 to be an axion and use only experimental data to determine its 
couplings. We point out that the self couplings are completely unconstrained 
experimentally. 

We propose the following mechanism for trapping pseudoscalars. We as¬ 
sume that the self coupling A is of order unity. The photons in the sun’s 
core convert into pseudoscalars due to the Primakoff process shown in fig. P 
The temperatures are such that almost the entire flux of photons emerging 
from the core is converted into pseudoscalars. The pseudoscalars begin ac¬ 
cumulating inside the core due to scattering on other pseudoscalars through 
the process 00 ^ 00. This is possible if the self coupling of pseudoscalars 
A is sufficiently strong. We can also have a signihcant cross section for the 
process 00 ^ 0000 through the loop diagram, £g. El if the coupling A is 
strong. This process leads to a degradation of the energy per particle for 
pseudoscalars. Hence as the pseudoscalars accumulate they start producing 
pseudoscalars of lower energy through this process and their density increases 
rapidly. This process stops when the energy per particle approaches the mass 
of the pseudoscalars. As we show the final density of pseudoscalars is suffi¬ 
cient to trap pseudoscalars inside the star. We also hnd that the mean free 
path of pseudoscalars becomes so small that radiative transport [B] occurs 
primarily through photons. 
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Figure 1: Feynman diagram for the Primakoff process for the conversion of 
photons 7 (fc) into pseudoscalars 



Figure 2: The process (jxj) —*• through a loop diagram. 



Figure 3: Feynman diagram for the inverse Primakoff process for the conver¬ 
sion of pseudoscalars (j){k) to photons 7(/c')- 
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We first compute the pseudoscalar density in the core of the sun assuming 
the current equilibrium conditions. We assume the core radius to be one-fifth 
of the solar radius Rq. The mean density in this region is approximately 
60 gm/cm^. The average temperature in the sun’s core is approximately 
1.7 X 10^ K. The pseudoscalars are produced by the Primakoff like process 
shown in fig. Q] and have mean energy roughly equal to 800 eV. We assume 
that the production rate of photons inside the core is equal to their rate 
of conversion into pseudoscalars. This is a reasonable assumption since the 
energy loss due to pseudoscalar production by this process is found to be 
larger than luminosity of sun. The solar luminosity is 3.9 x 10^^ erg/s. Hence 
the mean photon energy production rate inside the sun is roughly 2 erg/ (gm 
sec). It is clear that the production rate will be larger inside the core due 
to its higher temperature. For a conservative estimate we take the photon 
energy production rate inside the core to be 1 erg/(gm sec). The mean 
thermal energy per photon inside the core is roughly 800 eV 1.3 x 10“® 
ergs). Hence the number of photons produced per second is approximately 
5.21 X10^^. Since we have assumed that all the photons produced convert into 
pseudoscalars we find that this is also the number of pseudoscalars produced 
per second. Assuming that these particles escape freely from the sun we find 
that the flux of pseudoscalars from the sun’s core is 2.14 x lO^'^cm^^s”^. If 
the pseudoscalars are in a steady state the number density of pseudoscalars 
inside the core is 7.1 x 10®cm“^. 

We next determine the mean free path of a pseudoscalar propagating 
through this medium. The cross section of the process 0 (pi) 0 (p 2 ) —^ 0 (P 3 ) 0 (P 4 ) 
due to the (j)'^ coupling at leading order in perturbation theory is given by 




^ (fxf) 


327rE2 


= 9.94 X lO^A^ 


IQ-^ GeV 

Erm. 


GeV 


-2 


(3) 


With A = 1 and Ecm = 1 KeV, the mean free path = l/(n 0 (T<^ 0 ) = 
3.6 X 10^ cm. This is almost three orders of magnitude smaller than the 
core radius O. 2 i ?0 = 1.4 x 10^° cm. This implies that due to self interaction 
pseudoscalars will not be able to escape freely from the sun but will instead 
start accumulating inside the core. Hence our initial assumption that the 
pseudoscalars escape freely from the sun is not valid. 

We next consider the contribution due to higher order processes. We 
focus on the process 00 —> 0000 for which the leading order diagram is 
shown in £g. |21 This process is suppressed compared to the leading order 
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process by two powers of ax = A^/47r. With A = 1, this process is down 
only by a factor of 100 and can contribute signihcantly. This process leads to 
fragmentation of (j) particles and hence will degrade their energy per particle. 
The fragmentation stops when the energy per particle approaches their mass, 
which is of order 1 meV. Since the initial energy is 1 KeV, this will lead to an 
increase in the number of particles by a factor of 10® and hence will result in 
a considerable enhancement of the number density of (j) particles. The center 
of mass energy is now equal to x/2m^E^. This further enhances the cross 
section for the process 00 ^ 00 which becomes 2.5 x 10^® GeV“^. 

The above analysis shows that, during the main sequence phase, the 
pseudoscalars will accumulate inside the sun due to self interaction. As they 
accumulate, processes which deplete pseudoscalars also start contributing 
signihcantly. The dominant process contributing to depletion is the inverse 
Primakoff process shown in £g. El The mean energy of the pseudoscalars 
is much smaller than that of photons. Hence the photon produced by the 
inverse Primakoff process will have much lower energy, of the order of meV. 
The photons produced will be quickly absorbed by electrons by the inverse 
Bremsstrahlung process and thermalize in the medium. Eventually the en¬ 
ergy deposited into pseudoscalars will equal the energy released by pseu¬ 
doscalars due to their conversion back into photons. We next determine the 
density prohle of the pseudoscalars in this equilibrium situation. 

The cross section for the inverse process (hg. Ej) is approximately 10“^^ 
GeV“^ for « nie, assuming the target particles to be electrons. Their 
number density will reach steady state only when the rate at which energy 
is gained by the pseudoscalars becomes equal to the rate at which they loose 
energy. The mean energy of photons is about 10® times larger than that of 
pseudoscalars. Hence a pseudoscalar produced by the primakoff like process 
(hg. ID) will have roughly 10® times the energy of the photon produced by 
the inverse process (hg. El)- The production rate of pseudoscalars per unit 
volume is a^x^fpxn^rixv where n = c is the relative velocity and is the 
photon density of the medium. We set this equal to 10“® times the conversion 
rate of pseudoscalars into photons o'cpx^'yxnxn^v. This gives us the hnal 
number density of pseudoscalars when the system reaches steady state = 
10^n^a-yx^<j,x / o'tpx^-yx- The number density of photons in the core is roughly 
10^^ per cm^. Since the two cross sections are approximately equal, we hnd 
that the pseudoscalar number density required to achieve steady state is 
roughly 10^® per cm^. This leads to a mean free path of pseudoscalars of 
order 10“^^ cm inside the core of the sun. This is much smaller than the 
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mean free path of photons inside the core. Hence the contribution of 0 
particles to radiative transport inside sun would be negligible compared to 
the contribution of photons, once the sun becomes a main sequence star. 

So far we have only considered pseudoscalars inside the sun’s core as¬ 
suming steady state conditions. We expect the pseudoscalar density to be 
signihcant outside the core, extending much beyond the solar radius. This 
pseudoscalar halo will be gravitationally bound to the sun. We assume it to 
be spherically symmetric since we expect its rotational speed to be small. Due 
to the fragmentation process the pseudoscalar kinetic energies after emerging 
from the core can be atmost as large as their mass. Hence their velocity may 
still be comparable to the velocity of light, which is much larger in compari¬ 
son to the escape velocity from the surface of the sun, Vq ^ 6 x 10^ m/s. We 
assume that the pseudoscalars loose energy as they propagate inside the sun 
such that their velocities become nonrelativistic. There are many processes 
by which they can loose energy. For this we assume additional terms in the 
interaction lagrangian 


A-C/ = 5-^00750 + Ai0^ . (4) 

where the term proportional to Ai is parity violating. We are allowed to 
add such terms since parity is not a symmetry of nature. The coupling 
is also constrained by astrophysical processes. However, due to pseudoscalar 
self coupling and accumulation inside the sun the limits on are also much 
less stringent in comparison to the limit, g^ < 0.80 x 10“^° ITU], imposed by 
assuming that the pseudoscalar is an invisible axion. Laboratory bounds on 
this coupling also exist if we assume a very small mass of the pseudoscalar, 
i.e. if < 10“® eV HH. The coupling Ai is unconstrained. Some of 
the processes which can contribute to energy loss are shown in fig. 0 The 
diagram, fig. |3^, by itself can contribute significantly to energy loss if the 
coupling gA 10 which is not ruled out by any laboratory experiments. 
Furthermore the diagram, fig. EId, is only first order in the coupling gA and 
may lead to sizeable energy loss due to scattering on non-relativistic electrons. 

We may estimate the density prohle of the pseudoscalar halo for r > 
Rq by imposing steady state conditions. The pressure of pseudoscalars at 
distance r from center of sun is P = > /3. The mean velocity at 

distance r is given by 


2 2G{Mq + Mr) 

< >= C -I-^- - 

r 


(5) 
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Figure 4: Some of the diagrams which contribute to energy loss of pseu¬ 
doscalars as they propagate through sun. Here the pseudoscalar, electron 
and photon are represented by dashed, solid and wayy lines respectively. 
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where c =< Vq > —2G{Mq + Mr^)/Rq is a constant, Mq is the mass of 
the sun and Mr is the total mass of pseudoscalar particles within a sphere of 
radius r centered at the center of sun. We also have 


dMr 


dr 


and 


dP 

dr 


= ATTr^p^{r) , 


G{Mq + Mr)p 


Using equations El El and [7| and the equation P = p^ < > /3, we hnd, 


df^ 

dr 


<v^> 


G{Mq + Mr)p,j) 


+ SvrGpJr 


( 6 ) 

(7) 

( 8 ) 


Since the right hand side is negative for all r, we hnd that p^ decreases with 
r. In the small r regime where Mr « Mq and STip^j^r^ << Mq we hnd that 
p(j,{r) / P(j){ro) \/v^{r)/v^{ro). In the large r regime, Mr » Mq we hnd 
that oc 1/r^ and hence Mr oc r. 

In hg. Elwe show results of the numerical integration of coupled equations 
ElandlHl All masses are expressed in units of solar mass and distances in units 
of astronomical unit (AU). We have set the mass of the pseudoscalar equal 
to 1 meV. The results are shown for diherent values of the parameter c/G. 
We integrate the equations starting from the core of the sun. The equations 
given above can be easily generalized for Pq < r < i?core- For c/G > 0 
the pseudoscalars have enough energy to escape the gravitational attraction 
of the sun and hence the halo extends to inhnite distance. For values of 
c/G < —1 we hnd that the radius of the pseudoscalar halo is quite small and 
contributes negligibly to the mass of the solar system. Hence we hnd that 
for a wide range of parameters we can evade all the astrophysical bounds on 
light pseudoscalars. 

For c/G > 0, we hnd that Mr » Mq for r » 1 AU. Very large values 
of Mr are clearly unphysical and suggest that in this case the pseudoscalars 
never reach equilibrium. In such cases pseudoscalars may rapidly cool the 
star which may never reach the main sequence. We also hnd that the large 
r behavior, Mr » Mq, of the density of pseudoscalars is precisely the same 
as that required for the galactic dark matter. This raises the possibility that 
the galactic dark matter might be identihed with these light pseudoscalars. 
If these pseudoscalars populate the galaxy then they will also exert pressure 











Figure 5: The mass contribution of pseudoscalars, as a function of the 
distance from the center of the sun. 
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on the pseudoscalars emitted by the sun. This may also have considerable 
influence on the pseudoscalar density profile and will be investigated in a 
future publication. 

So far we have shown that, assuming the current equilibrium conditions 
inside sun’s core, the pseudoscalars contribute negligibly to radiative trans¬ 
port. Since all the bounds imposed in the literature consider perturbations 
around the equilibrium conditions, our analysis show that such bounds are 
not valid in the presence of self couplings of pseudoscalars. 

We next consider the approach of the star to the main sequence phase. 
Since we are considering a relatively large pseudoscalar photon coupling the 
pseudoscalar flux can be significant even when the star has not reached the 
main sequence. When the temperature is lower than the current temperature 
the pseudoscalar flux and hence their density inside the core will be smaller. 
Hence for a certain temperature range the pseudoscalars will escape freely. 
We consider the phase when the star undergoes quasi-static collapse due to its 
gravitational attraction. By virial theorem the total energy of a star in equi¬ 
librium is equal to half its potential energy. Hence its energy loss is limited 
by virial theorem and pseudoscalar production cannot change its evolution 
drastically. We point out that the star will evolve quasi-statically even in the 
presence of pseudoscalars since the processes leading to equilibrium proceed 
at a much faster rate in comparison to the production rate of pseudoscalars. 
As the core temperature increases, the rate of photon and hence pseudoscalar 
production also increases. Once the photon production rate reaches values 
close to what is observed today, pseudoscalars will start accumulating and 
their contribution to radiative transport will rapidly decrease. 


3 Conclusions 

In conclusion we And that it is possible to evade all astrophysical limits on 
pseudoscalar photon coupling if the self coupling of pseudoscalars is suffi¬ 
ciently strong. The pseudoscalars start accumulating inside the sun due to 
scattering with other pseudoscalars. The energy per pseudoscalar degrades 
due to higher order processes as well as by energy loss due to scattering 
on electrons, nucleons and pseudoscalars. The pseudoscalar number density 
eventually reaches steady state due to conversion of pseudoscalars back into 
photons. We find that the mean free path of these particles inside the core 
is much smaller in comparison to the mean free path of photons. Hence 
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they contribute almost negligibly to radiative transport. In our analysis we 
have used the pseudoscalar coupling and mass extracted by the PVLAS ex¬ 
periment. Our results show that the results of the PVLAS experiment are 
consistent with astrophysical limits if we allow strong self couplings. 
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